Late-time decay of scalar perturbations outside rotating black holes 
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We present an analytic method for calculating the late-time tails of a linear scalar field outside a 
Kerr black hole. We give the asymptotic behavior at timelike infinity (for fixed r), at future null 
infinity, and along the event horizon (EH) . In all three asymptotic regions we find a power-law decay. 
We show that the power indices describing the decay of the various modes at fixed r differ from the 
corresponding Schwarzschild values. Also, the scalar field oscillates along the null generators of the 
EH (with advanced-time frequency proportional to the mode's magnetic number m). 
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The no hair principle for black holes implies that grav- 
itational field outside a generic black hole relaxes at late 
time to the stationary Kerr-Newman geometry. For lin- 
ear test fields (either scalar, electromagnetic, or gravi- 
tational) outside a spherically-symmetric Schwarzschild 
black hole, it was shown by Price JjJ that all radiative 
multipoles die off at late time with a t~ 2l ~ 3 power-law 
tail 0], where I is the mode's multipole number, and t 
is the Schwarzschild time coordinate. Later, this result 
was confirmed using several different techniques, both 
analytic and numerical &-k| . The relevance of the per- 
turbative (linear) results to the fully nonlinear late-time 
behavior was demonstrated in numerical simulations of 
a fully non- linear, self- gravitating, spherically-symmetric 
scalar field [@j|. 

It is well known, however, that realistic astrophysical 
black holes are spinning and not spherically-symmetric 
M. Therefore, for astrophysical applications it is ex- 
tremely important to generalize the above analyses from 
the Schwarzschild background to the more realistic Kerr 
background. A first progress in this direction has been 
achieved recently with the numerical simulation of linear 
fields in the Kerr background, by Krivan et al. M,0|. 
Yet, a thorough analytic treatment of this problem has 
not been carried out so far (l"^|l3| ]. 

The goal of this Letter is to present an analytic method 
for calculating the late time behavior of a linear mass- 
less scalar field outside a Kerr black hole. This method 
was recently applied to the simpler Schwarzschild case 
as a test-bed |6|, in which case the well known late-time 
inverse-power tails were reproduced. In this Letter we 
outline the application of this method to the Kerr case, 
and present the main results. In particular, we calcu- 
late the power indices characterizing the late-time de- 
cay of the various modes at future null infinity, at fixed 
r, and at the EH. Quite interestingly, we find that at 
fixed r these indices are different than those found in the 
Schwarzschild case. Full details of the calculations are 
given in Ref. jl5|]. Throughout this paper we shall use 
the standard Boyer-Lindquist coordinates t,r,9,tp, and 
relativistic units C = G = 1. 

The main difficulty in analyzing perturbations over a 



Kerr background is the nontrivial dependence on 9. The 
separation of variables by the Teukolsky equation pl| is 
only applicable to the Fourier-decomposed field, because 
the spheroidal harmonics used for the separation of the 9, 
ip variables explicitly depend on the temporal frequency 
u>. The final goal is to calculate the late-time decay of the 
field, along with its angular dependence, in terms of the 
time t. Obviously, an expression of this angular depen- 
dence in terms of the (w-dependent) spheroidal harmon- 
ics would be useless. This motivates one to carry out the 
entire analysis in terms of the spherical harmonics iy™. 
The difficulty is, however, that due to the breakdown of 
spherical symmetry, modes of different I (but the same 
m) are coupled; Namely, there are "interactions" between 
modes. The main challenge is to handle this interaction 
and to analyze its effect on the late-time decay. 

In principle, it is possible to carry out the analysis in 
the frequency domain, and then Fourier-integrate over 
all frequencies to recover the late-time behavior in the 
time domain, as was done in the Schwarzschild case [0,9 . 
We find it very difficult, however, to properly apply this 
method to the present problem, due to the following rea- 
son. In the Schwarzschild case the analysis can be much 
simplified by taking the limit uj — > 0. In the Kerr case, 
this limit docs not lead to a correct description of the 
interaction between modes |13|. As we show below, it is 
this interaction which dominates the late-time decay of 
the modes I > 2 at fixed r. In the limit uj — > one simply 
(and incorrectly) recovers the Schwarzschild result |p7| . 

In view of these considerations, we found it simpler 
to carry out the entire analysis in the time domain (i.e. 
without a Fourier decomposition). To overcome the dif- 
ficulties caused by the interaction between modes, in the 
first part of the analysis we use an iteration scheme in 
which we iterate over the interaction term (along with 
the other curvature-induced terms in the field equation). 
In the second part of the analysis we use the late-time 
expansion, which is essentially an expansion in inverse 
powers of advanced time. (Both methods are generaliza- 
tions of those used in Refs. UqM for the spherical case). 

The Klein- Gordon field $ in Kerr geometry obeys 
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where M and a are, correspondingly, the black hole's 
mass and specific angular momentum, and A = r 2 — 
2Mr + a 2 . Decomposing $ into spherical harmonics in 
the standard manner, $ = J2i m 4> lm {t, r ) ^ m (#, <f), an d 
defining w m = y/r 2 + a 2 <p lm , the original field equation 
(fy) becomes (for each m) 
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where cq, c_ and c+ are certain coefficients depending on 
I and m (w: 
and where 



I and m (with c m ' = 0; no other coefficients vanish), 
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The coordinates u and u are defined by u = t — r* and 
w = £ + r 1 *, with r*(r) obeying dr*/dr = (a 2 + r 2 )/A. 
[In Eq. (0), and also in most of the equations below, we 
omit the index m for brevity. Note that due to the axial 
symmetry, modes with different m do not interact.] 

Equation (0) is an infinite set of coupled equations for 
the various modes of the field, with the last two terms in 
the square brackets describing the interactions between 
modes of different I. 

The set-up of initial data for the evolution problem 
is similar to that used in Ref. [0 for the Schwarzschild 
case [see figure 2 and Eq. (7) therein]. That is, we as- 
sume that $ is specified along a pair of hypersurfaces 
v = vq and u — uq, and without loss of generality we 
take vq = 0. For convenience we consider a situation of 
an outgoing pulse at v = 0, which starts immediately 
after the outgoing ray u = Uq. Wc further assume that 
— Uq 3> M, and that the pulse is arbitrarily-shaped but 
relatively short, which considerably simplifies the analy- 
sis [|18lJ6|]. Although this type of initial data is not the 
most general one, it is nevertheless reasonably generic, 
and we expect the resultant asymptotic behavior to be 
characteristic of the generic situation. We also assume 
here that the initial outgoing pulse has a rather generic 
angular distribution, so it includes all the spherical har- 
monics (and especially the component 1 = 0). (Later we 
shall briefly discuss the more subtle case, in which no 
I = component is present at the initial data.) 

To evolve these initial data and analyze the late-time 
behavior, we shall proceed in two steps. In the first step, 
we calculate the late time (i.e. u ^> M) form of the field 
at future null-infinity (v —> oo). In the second step we 
derive an expression for the field at a fixed r at t > M 
(and also along the EH at n> M). 



Asymptotic behavior at future null- infinity. We now 
apply the iteration scheme introduced in Refs. Eq,0|, and 
decompose ^ lm as 
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The components "J^v* arc defined by the hierarchy of 
equations 
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where S l = and (for N > 1) 
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*' and ^ l N>1 = 



along with the initial conditions \& 
at v = and u = Uq. Here, Vg(r*) is the Minkowski- like 
potential defined in Ref. || as a function of r* [Eqs. (8,60) 
therein] , and SV^ (r) = V^ — Vq . Formal summation over 
N recovers the original field equation and initial data for 
the complete fields if? 1 . 

The field equation (g), together with the above initial 
conditions, constitutes a hierarchy of initial- value prob- 
lems for the various functions if? 1 ^ 1 , which, in principle, 
we may solve one by one (first for N = 0, then for N = 1, 
etc.). Notice that the potential Vb(r„) (and hence the en- 
tire N = equation) is independent of the spin parame- 
ter a. The solution of this equation, the function if? l , is 
given explicitly in Ref. B (see section IV therein). This 
function decays exponentially at late time, so it does not 
contribute to the power-law tail. Rather, it serves as a 
source for terms ^jv>i, which do provide power-law tails. 
For each N > 1, the field equation can be formally solved 
in terms of a Green's function: 
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where G (u, v;u',v') is the (retarded) Green's function 
associated with the zero-order operator d u d v + Vq. An 
analytic expression for G was derived in section V of Ref. 
H . This, in principle, enables the solution of the field 
equation (||) for all N and I. 

The functions if?\ will primarily concern us here, be- 
cause it is the term N — 1 which dominates the late-time 
behavior of if? at null infinity in the generic situation. 
It is convenient to consider separately the contributions 
from the various terms in Eq. (0) [through Eq. (pi)] to 
if?\. Consider first the contribution from the term pro- 
portional to SVft. This potential can be expressed as 
b~V l K = SVg + 5V* , where 6Vg is the (a-independent) cor- 
responding Schwarzschild contribution, and SV^ is an a- 
dependent correction term. A direct calculation shows 



that at r > M, 5V* decays faster than SVg by a fac- 
tor proportional to a 2 /(Mr*). An explicit evaluation of 
the integral in Eq. (J^) then yields that this extra factor 
leads to an extra u _1 factor in the late-time asymptotic 
behavior of ^[ at null infinity (L5J. Thus, the dominant 
contribution to W^ from SV^ at null infinity, which we 
denote by W 1: is the same as in the Schwarzschild case 
(cf. Eq. (58) inRef. §): 
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where Ai is given explicitly in H as a linear functional 
of the Z-component of the initial pulse . 

Consider next the contribution to ^[ due to the other 
part of S{, i.e. the part containing ^-derivatives in Eq. 
(pf) . A direct evaluation of the integrals in Eq. (R) shows 
that the late-time contribution of this part at null infin- 
ity is proportional to u~ l ~ 3 or smaller |15), and is hence 
negligible. The only exception is the contribution from 
^ott (for I > 2), which is proportional to u~ l ~ 2 too, so 
it does not cause a qualitative change in the asymptotic 
decay (||). Moreover, the coefficient of this term is re- 
duced by a factor proportional to (a/u ) 2 <C 1 compared 
to A{ , so the overall tail amplitudes are still given by Eq. 
(||) at the leading order. 

We still need to consider the terms N > 2. A com- 
plete analysis of these terms has not been carried out 
yet. In the Schwarzschild case, simple considerations sug- 
gest that at null infinity all these terms decay like u~ l ~ 2 , 
though with coefficients smaller than that of \&i by a 
factor (M/wo)^^ 1 . (This was also verified by numeri- 
cal simulations Q, which also suggested convergence of 
the sum (H) at null infinity.) Hence the N > 2 terms 
do not alter the power index, and, moreover, in the case 
— uq 3> M considered here, they do not significantly af- 
fect the coefficients. All these considerations apply to 
the Kerr case as well p5|,E9l. Assuming that the terms 
N > 2 indeed behave in that manner, we find that at 
late time, the scalar field at null infinity is given by 
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(null infinity, u 3> M), 



(9) 



where the coefficients Ai coincide with those of Eq. (g) 
to the leading order in M/uq. 

Derivation of $ at r=const: the late-time expansion. 
We now derive an expression for the late-time behavior at 
any fixed r outside the black hole and along the EH, ac- 
curate to the leading order in M/t or M/v, respectively. 
To that end we employ the late-time expansion used in 
the Schwarzschild case (cf. Ref. [pi): 
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fc=0 



As it turns out, this expansion is consistent with the field 
equation, with the regularity condition at the EH, and 
with the form of the field at null-infinity The parameter 



fco > is by definition ^-independent, and will later be 
determined through matching to null infinity. (For I > 
some of the first terms in the sum (|I(j) vanish, as will 
become apparent below.) 

Substituting the expansion ( |lCj| ) in the field equation 
(0), and collecting terms of the same power in v, the par- 
tial differential equation becomes an infinite coupled set 
of ordinary equations for the unknown functions i^(r), 
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where a prime denotes d/dr (and where F l k , <0 = 0). 
Here, too, the source term Z l k exhibits interactions with 
modes V ^ I. However, since Z l k only depends on terms 
Fj., with k' < k, the system (|ll| ) is effectively decoupled, 
as the various ordinary equations may be solved one at 
a time. It is possible to formally write down the general 
solution for any F k via the Green's- function method [[15| . 
Consider first the term k = 0, which dominates the over- 
all late-time asymptotic behavior at fixed r. The function 
-F|. =0 satisfies a homogeneous equation (actually the sta- 
tionary field equation), whose general solution is given 
by F l Q = aiP~~<(p) + 6 ; P / +7 ( / o). Here, a/ and h are (yet) 
arbitrary constants, p = (2r — r + — r—)/(r+ — r_) where 
r± = M ± (M 2 — a 2 ) , and P l 1 are the two complex- 
conjugated associated Legendre functions of the first kind 
p0| with an imaginary index 7 = im [2a/ (r+ — r_)]. The 
functions P l 7 (which are special cases of the Hypergeo- 
metric function) have the form 

P^{P) = Vf{p) x [(p + \)/{p I )] ±7/2 , (12) 

in which V t 7 are (complex) polynomials of order I (non- 
vanishing at r — * r+) . For m j^ these functions oscillate 
rapidly towards the EH (r -^ r +1 p — ► 1), 

P z ±7 (r -> r+) oc (p - 1) T7/2 oc exp(=Rmn + r*), (13) 

where Q, + = a/(2Mr+). 

One of the two coefficients a; , 6; is to be determined 
from the regularity condition at the EH. Here one must 
recall that the Boyer-Lindquist coordinate tp is singular 
at the EH. Using the regularized azimuthal coordinate 
ip + = ip — il + t pi[ instead, one finds 
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Since the factor in square brackets is regular at the EH 
(but the next factor is not), it follows from the regularity 
condition that bi = 0, hence Fg = aiP t r (p). 

The parameter a; is to be determined through match- 
ing to null infinity. The asymptotic form of Fq as r — > 00 
is Fy(r ;§> M) oc r l . Substitution in Eq. ( |l"o| ) (taking into 
account the contribution of the terms k > as well, which 



are not negligible at null infinity, as explained in || ) , one 
obtains at null infinity ^f l ex aiu 



l+l-ko 



m 



Matching 

this expression to Eq. (||) for I = yields ko = 3. This 
value of fco yields a consistent matching for any I, im- 
plying a;>i = (that is, the modes I > 1 are excited 
only at k > 0). One finds that the dominant mode I = 
decays like v~ 3 ex £~ 3 at fixed r (and large t), as in the 
Schwarzschild case. 

The interaction between modes has a crucial effect on 
the decay rate of modes I > 2. Without this interaction, 
one can verify that a mode I, m is excited at k = 21, lead- 
ing to a decay rate t~ 2l ~ 3 (as in the Schwarzschild case) 
p2| . The interaction changes this situation. Consider, 
for example, the mode I = 2, m = 0. This mode has 
a vanishing source term Z]^[ (as F^z^ = 0), and one 
can show (using the argument of [p2[ ) that Ff.z^(r) = 0. 
On the other hand, at k = 2 there is a non- vanishing 
source term Zjz^ ex C-Fj^, which necessarily leads to 
a non- vanishing function Fjzz^r). Thus, the decay rate 
of this mode at fixed r is v~ k °~ k = v~ 5 ex £~ 5 , which 
differs from the corresponding Schwarzschild rate, t . 
This simple consideration is easily extended to all other 
modes m,l, and one finds |l5| 
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(Xt-i-M-3-q ( nxec i rj t>M,|r*|), (15) 



where q = for even I + m and q — 1 for odd I + m. The 
late-time behavior of a mode l,m at the EH (expressed 
in regular coordinates) is found to be 

¥ m Y? n {6, ip) ex y ; m (6», (p+) e lmQ + v tj-'-NI-s-s. (i 6 ) 

(This power index may be changed if the relevant func- 
tion Fj,(r) happens to vanish as r — ► r + .) 

In summary, the late-time behavior of the various 
modes in the three asymptotic regions is given in Eqs. 
(i)> (111), and ©• (For the dominant mode Z = 0, the 
amplitude coefficients in all three asymptotic regions are 
given explicitly in Ref. |l5J].) Our analysis indicates two 
phenomena special to the Kerr case: 

A. Oscillations along the EH — cf. Eq. ( pL6|) . 

B. The interaction between modes: Due to this interac- 
tion, the power index at fixed r is I + \m\ + 3 + q. (When 
the interactions are ignored, one obtains the standard 
Schwarzschild exponent 21 + 3 [ft7|.) 

Throughout this paper we have assumed that the ini- 
tial pulse includes all the modes (and in particular, the 
dominant mode 1 = 0). In non-generic situations in which 
the \ow-l modes are absent at the initial data, the interac- 
tion between modes may dominate the overall late-time 
behavior already at null infinity. For example, assume 
that the initial pulse is a pure I — 2, m = mode. Then, 
without the interactions, at null infinity $ would be dom- 
inated by 'J^! oc u~ A . However, the interaction excites 
(at N — 2) an I — mode with a u~ 2 tail, that dominates 
the late-time behavior. This behavior has been observed 



numerically by Krivan et al. [M. In Ref. |L5| this phe- 
nomenon will be discussed in more detail, along with its 
implications to the asymptotic behavior at fixed r. 

We should emphasize that despite the relative simplic- 
ity of the calculation scheme presented here, the mathe- 
matical question of convergence of the various expansions 
involved is still open (though there is evidence for con- 
vergence) . This is the situation even in the Schwarzschild 
case. Additional mathematical subtleties arise in the 
Kerr case, which we further discuss in p5[ . 

We finally note that numerical simulations pi are con- 
sistent with our analytic results for the power indices of 
the overall perturbation at fixed r |23] . It will be interest- 
ing to numerically test our prediction of the power index 
I + \m\ + 3 + q for the individual I, m modes at fixed r. 

Note added: After this manuscript has been submit- 
ted, Hod analyzed the mode coupling in Kerr spacetime 
|L4j . In the scalar field case he recovers our results. 
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